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The effect of Coulomb and short-range interactions on the spectral properties of two-dimensional 
disordered systems with two spinless fermions is investigated by numerical scaling techniques. The 
size independent universality of the critical nearest level-spacing distribution P(s) allows one to 
find a delocalization transition at a critical disorder W c for any non-zero value of the interaction 
strength. At the critical point the spacings distribution has a small- s behavior P c (s) oc s, and a 
Poisson-like decay at large spacings. 

PACS number(s): 71.30, 72.15 Rn, 71.55 Jv 

It has been shown that the statistical properties of 
spectra of disordered one-electron systems are closely re- 
lated to the localization properties of the corresponding 
wave functions In the localized regime, states close 

in energy have an exponentially small overlap and their 
levels are uncorrelated. So the corresponding normalized 
spacings s are distributed according to the Poisson law 



Pp(s) = exp(-s) 



(1) 



On the other hand, in the metallic regime the large 
overlap of delocalized states induces correlations in the 
spectrum leading to the well known level repulsion ef- 
fect. If the system is invariant under rotations and un- 
der time-reversal symmetry the spacings follow Wigner- 
Dyson statistics 



7T / 7T , 

Pw{s) = 2 Sex P V ~4 S ' 



(2) 



The two limiting behaviors given by Eqs. (pj) and (^J) cor- 
respond to an infinite system. In Refs. |2|.[1| was found 
that a third, size independent distribution, appears at 
the metal-insulator transition (MIT) in the three dimen- 
sional (3D) Anderson model. This behavior is the basis 
of a different method to detect the MIT whose advan- 
tage resides in that no knowledge of the eigenfunctions 
is needed. We will extend this procedure to two-electron 
interacting disordered systems. 

The study of the energy spectrum has already pro- 
vided important results about the nature of the states 
and about the transport regimes that one can expect. 
In fact, the method has been successful in the location of 
the MIT in the 3D Anderson model > m two dimen- 

sional (2D) disordered systems with symplectic symme- 
try [p|-pl| , in the determination of the critical interaction 
strength in many-body systems |l2|-|l4j , and in the study 
of transport regimes in the 2D Anderson model Ollq ]. 

As dimension two is marginal, one can expect a 
interaction-driven delocalization transition in a 2D in- 
teracting disordered system because of the following rea- 
sons. First of all, the scaling theory of localization includ- 
ing the combined effects of disorder and interactions pre- 
dicts that a 2D system may remain metallic even in the 



limit of zero temperature |17 lSy . Secondly, many recent 
experimental results have presented strong evidence for 
a delocalization transition in 2D systems without mag- 
netic scattering effects employing different materials and 
designs |l9| -^5| . All these experiments show clear indica- 
tions that strong electron-electron interactions partially 
suppress the quantum interference effects responsible for 
localization. Furthermore, a recent numerical study of 
two-interacting bosons in 2D by means of their decay 
length in long bars [^6| showed evidence for a localized 
to extended states transition. 

In this Letter we study the level statistics of two inter- 
acting electrons in a 2D random potential. We consider 
both long-range Coulomb interactions and short-range 
interactions and calculate numerically the correlations in 
the exact two-electron spectrum. Our main result is that 
the nearest level-spacing distribution P(s) is universal in 
the critical region showing a delocalization transition for 
any finite value of the interaction strength. 

We consider spinless fermions on a sample of size Lx L 
described by the standard Anderson-Hubbard Hamilto- 
nian 

H = t^2(4,j+l a i,j + 4+l,j a i,3 + h - c -) 

+ ^2 e i,j a i,j a i,j + H int , (3) 



where the operator a\ - (ciij) creates (destroys) an elec- 
tron at site of a square lattice and tij is the en- 
ergy of this site chosen randomly between (— W/2, W/2) 
with uniform probability. The hopping matrix element t 
is taken equal to —1 and the lattice constant equal to 1, 
which sets the energy and length scales, respectively. For 
the long-range case, the interaction Hamiltonian is given 



U 



i,j>k,l 



4,j a i,3 a k,l a k,l 
J- r k,l\ 



(4) 



while for the short-range case we choose a nearest- 
neighbor interaction Hamiltonian 



1 



Hint = U ^2(4,j<Hj4,j+l a ij+l 
+ a i,3 a iJ a i+l,3 a i+l,j) ■ 



V(L,W) 



(5) 



In both cases, in order to reduce edge effects, we use pe- 
riodic boundary conditions, i.e., 4+1 j ~ 4j' 4 l+i = 

4i' a lj =4,ji 4,0 =4,L- 

The L 4 eigenstates of the Hamiltonian (||) are either 
symmetric or antisymmetric with respect to the inter- 
change of the electron positions. We restrict our inves- 
tigation to the two-electron Hilbert subspace spanned in 
the basis of N = L 2 (L 2 — l)/2 antisymmetric products 
of one-electron states 



\4>i,j;k,l) 



1 



(4,3 a U 



z f at 

L k,l a i,3 



)|0) 



(6) 



where |0) is the vacuum state. The resulting N x N 
Hamiltonian matrix is numerically diagonalized using 
techniques for large sparse matrices, in particular a 
Lanczos tridiagonalization without reorthogonalization 
method |27| . The strength U of the interactions is varied 
between 0—10 and the disorder strength ranges the in- 
terval (0.2, 16). The system size varies from L = 6 to 20 
and the number of random realizations is such that for a 
given triad of {U, L, W} the number of studied eigenval- 
ues in the range (—1, 1) was kept around 2.5 x 10 4 . 

In order to minimize fluctuations in P(s), we used in 
the calculations for U > 1 the levels contained in the 
interval (—1,1) around the band center. Although this 
energy window exceeds the one-electron level spacing Ai, 
does not introduce hidden correlations (unlike for the 
non-interacting case U = or for U <C 1) in the two- 
electron spectrum because of the strong mixing of the 
states directly coupled by the interaction. In Ref. |2£| was 
shown, modeling the full Hamiltonian (^) by a Gaussian 
ensemble of random matrices with a preferential basis 
p9[ , p0| , that the spread width T 2 of a two-electron state 
is of the order of A^ In any case, we verified that, within 
the statistical errors, there is no significant changes in 
P(s) if we used an energy window A e < Ai. 

To quantitatively analyze the scaling of the statisti- 
cal properties of the eigenvalue spectrum one can use 
one of the two following variables widely employed. 
The first one is given by j) = f™[P(s) - 

Pw(s)]ds/ f s [Pp(s) — Pw(s)]ds, where for sq one can 
choose one of the two crossing points of the curves Pp{s) 
and Pw(s), namely s Q = 0.473 and s = 2.002. The for- 
mer is representative of the small- s part of P(s), whereas 
the latter is representative for its tail. The second vari- 
able |||y[l[ , based in the A3 statistics, which take into 
account correlations between non adjacent levels, is given 
V = M Jo 30 A.a(L)dL. In this Letter we use a different 
scaling variable [[l6|,[32j based on the variance of P(s), 
since it contains information about the whole distribu- 
tion, namely 



var(s) - 0.273 
1 - 0.273 



(7) 



which describes the relative deviation of var(s) from the 
Wigner-Dyson limit due to the finiteness of the system. 
In Eq. (@) var(s) = (s 2 ) - (s) 2 , and 0.273 and 1 are 
the variances of Wigner-Dyson and Poisson distributions, 
respectively. In this way, r](L, W) ranges from [P{s) = 
Pw] to 1 [P(s) = Pp]. We have checked that our results 
are basically independent of which of the three definitions 
of 77 mentioned is used, although the results based on Eq. 
(0) are more accurate. 

Figure 1 shows the disorder dependence of 77 for a long- 
range interaction with U = 1 for different system sizes: 
L = 6 (o), 9 (□), 16 (o) and 20 (y). One can distinguish 
three different regions, in which r\ increases, remain con- 
stant, or decreases with increasing L. The curves for 
different sizes cross at a common point, W c , which cor- 
responds to the interaction-induced derealization tran- 
sition. 




16 



FIG. 1. Scaling variable 77 as a function of W for a 
long-range interaction with U = 1. System sizes: L = 6 (o), 
9 (□), 16 (o) and 20 (v)- The critical disorder is given by 
the point for which r\ is independent of L. Upper inset: Veri- 
fication of the one-parameter scaling assumption for r\{L, W). 
The horizontal scale is arbitrary. Lower inset: U dependence 

Of We- 
ill the following we check whether r/(L, W) can be ex- 
pressed by a one-parameter scaling function 



n(L,w) = f(L/aw)) , 



(8) 



where the scaling parameter £ is the two-electron localiza- 
tion length in the localized regime, and the two-particle 
correlation length in the extended regime. To perform 
this scaling procedure with the r\ data, the range of rj 
values for different L at any given disorder W\ must over- 
lap the range of 77 values for various L for at least one 
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different disorder W2 ■ Equation (|J) implies that in a log- 
log plot of r)(L, W) versus L all data should collapse in 
a common curve when translated by an amount ln£(VF) 
along the horizontal axis. This curve has a single branch 
when there is no transition, while it develops two sepa- 
rate branches when a transition is present. In the upper 
inset of Fig. 1, we re-plot the data of Fig. 1 by rescaling L. 
Note that one can fit most of the points onto a common 
curve with two branches, one growing (W > W c ) and 
another decaying (W < W c ). This results confirm the 
scaling hypothesis for r](L, W) and thus we have found 
the function £(W) up to an arbitrary factor. 

Next we determine the critical disorder W c and the 
critical exponent v. For this purpose we used the singu- 
larity of the correlation length near W c 



^{W) = io\W -W c \-» 



(9) 



where £0 is a constant. Using Eq. (||) and taking into 
account that W) is analytical for a finite system, we 
can write around the critical point 



(10) 



from which one can extract W c and v. In practice, we 
have truncated the series at n = 4. We performed an 
statistical analysis of the data in the range 10 < W < 
14 with the Levenberg-Marquardt method for nonlinear 
least-squares models. The most likely fit is determined 



by minimizing the \ 2 statistic of the fitting function (10). 
The critical disorder found is W c = 11.8 ± 0.2, and the 
corresponding critical exponent is equal to v = 1.2 ± 0.2. 

We repeated the same kind of calculations for different 
values of the interaction strength U including the non- 
interacting case (U = 0). For U <C 1 we obtained 77 using 
Eq. (0) m a energy window A e < Ai. The U dependence 
of the critical disorder W c is shown in the lower inset of 
Fig. 1. The W c found increases monotonously with U 
starting from very small values. In the non- interacting 
limit there is no critical behavior and all states are local- 
ized, as expected. 

We have also analyzed the behavior of two electrons 
subjected to a short-range interaction, Eq. (^). The cor- 
responding results are depicted if Fig. 2. The system 
sizes used are the same as those appearing in Fig. 1. 
For U — 1 the critical disorder and exponent found are 
W c = 10.4 ± 0.2 and 1.3 ± 0.2, respectively. The de- 
pendence of W c on U, shown in the lower inset of Fig. 
2, is slightly different than that found for the long-range 
case. W c increases with U and reaches a saturation value 
around U ~ 2. After saturation W c decreases for the 
largest value of the interaction strength used, U = 10. 
The reason is that for this value of U, as we have checked, 
the interaction slightly splits the two-interacting particle 
band into lower and upper Hubbard bands. 




FIG. 2. Same as in Fig. 1 for a short-range interaction with 
17 = 1. 

Finally we focus on the level statistics at the criti- 
cal point. For the orthogonal symmetry P c (s) oc s as 
s — > J2|,|33| , while for s 1 , two different analytical ex- 




assumes that 
has the asymp- 



pression were proposed. One of them 
P c (s) is Poissonian and the other 
totic form P c (s) oc exp(— A c s a ), where a is related to the 
critical exponent v and to the dimensionality d through 
a = 1 + [dv) . In order to diminish the magnitude 
of the relative fluctuations and to analyze the asymp- 
totic behavior in detail, it is more convenient to con- 
sider the cumulative level spacing distribution function 
I(s) = J s °° P(s')ds'. The Poisson distribution Eq. (@) 
and the Wigner surmise Eq. ^ yield Ip(s) — exp(— s) 
and Iw(s) = cxp(— 7rs 2 /4), respectively. 

We have checked that for various L considered at the 
critical disorder W c for a given U, the best fit using a x 2 
criterion in the interval 2 < s < 5 yields a w 1.0±0.1 and, 
therefore, P c (s) at large s is very close to a Poissonian de- 
cay, thus confirming the ideas of ||J4|]. A similar conclu- 
sion was obtained in Ref. Q in the Anderson transition 
using very large system sizes. The results for the nu- 
merical calculations of the critical I c (s) for a long-range 
interaction with U = 1 are shown in Fig. 3 for L = 12 
(o), 16 (□) and 20 (o). Note that the critical I c (s) is a 
L-independent universal scale-invariant function that in- 
terpolates between Wigner and Poisson limits. Solid and 
dashed lines correspond to Ip(s) and iw(s)j respectively. 
In the upper inset we display our results for I c (s) in the 
limit s ^> 1. One clearly see that the Poissonian tail is 
recovered for large spacings. This behavior is similar to 
that in the insulating regime, although the decay rate A c 
is larger than unity. The straight line fitting the data in 
the interval 2 < s < 5 is \nl c (s) = -1.32s + 0.33. In 
the lower inset we show the limiting behavior of P c (s) 
as s — > 0. Here we find that P c (s) oc s according with 
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predictions of Refs. |^|33|]. The straight line fitting the 
data in the interval 4 x 10~ 3 < s < 1CT 1 is P c (s) = 3s. 




FIG. 3. Integrated probability 7 c (s) at the critical disorder 
for a long-range interaction with {7 = 1. System size: L — 12 
(o), 16 (□), 20 (o). Solid and dashed lines are Ip(s) and 7w(s), 
respectively. Upper inset: Large-s part of 7 c (s); the fitted line 
is ln7 c (s) = —1.32s + 0.33. Lower inset: double logarithmic 
plot of P c (s) for s < 1; the fitted line is P c (s) = 3s. 

To summarize, we have investigated the energy levels 
statistics of a two-electron disordered 2D system with a 
long-range or a short-range interaction. In both cases, 
the scaling function was calculated using the statistics 
of the nearest-neighbor eigenvalues of the two-electron 
states, that have been obtained numerically. By per- 
forming a finite size scaling analysis we found a critical 
disorder W c , for any non-zero value of the interaction 
strength U. When 7/ — > 0, W c — > and so we recover 
the expected results of the non-interacting system where 
all the states are localized. P c (s) is linear for small s 
while the large-s part of P c (s) obtained is shown to have 
a Poisson-like decay. 
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